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ABSTRACT 

We report a study of the stellar content of the Near-infrared cluster [DBS2003] 157 embed- 
ded in the extended H ii region GAL 331.31-00.34, which is associated with the IRAS source 
16085-5138. JHK photometry was carried out in order to identify potential ionizing candi- 
dates, and the follow-up NIR spectroscopy allowed the spectral classification of some sources, 
including two O-type stars. A combination of NIR photometry and spectroscopy data was used 
to obtain the distance of these two stars, with the method of spectroscopic parallax: IRS 298 
(06 V, 3.35 + 0.61 kpc) and IRS 339 (09 V, 3.24 + 0.56 kpc). Adopting the average distance 
of 3.29 + 0.58 kpc and comparing the Lyman continuum luminosity of these stars with that 
required to account for the radio continuum flux of the H ii region, we conclude that these 
two stars are the ionizing sources of GAL331. 31-00. 34. Young stellar objects (YSOs) were 
searched by using our NIR photometry and MIR data from the GLIMPSE survey. The analysis 
of NIR and MIR colour-colour diagrams resulted in 47 YSO candidates. 

The GLIMPSE counterpart of IRAS 16085-5138, which presents IRAS colour indices 
compatible with an ultra-compact H ii region, has been identified. The analysis of its spectral 
energy distribution between 2 and 100/im revealed that this source shows a spectral index 
a = 3.6 between 2 and 25 /jm, which is typical of a YSO immersed in a protostellar envelope. 
Lower limits to the bolometric luminosity and the mass of the embedded protostar have been 
estimated as L = 7.7 x 1O^L0 and M = 10 Mq, respectively, which corresponds to a BO-B 1 V 
ZAMS star 

Key words: stars: early-type - H n regions - stars: pre-main-sequence. 



1 INTRODUCTION 

The study of the stellar content and the determination of the 
distances of Hii regions and star-forming complexes associated 
with massive molecular clouds are fundamental for the determi- 
nation of the spiral structure and the rotation curve of the Galaxy. 
Also, the assessment of Galactic gradients of chemical abundances 
and electron temperatures have strong dependence on the accu- 
racy of the distance estimates. The Norma region is a very inter- 
esting sector of the Galaxy for this investigation, since the line 



* Based on observations carried at the SOAR observatory, a joint project of 
the Ministerio de Ciencia, Tecnologia e Inovafao (MCTI) of the Republica 
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(NOAO), the University of North Carolina at Chapel Hill (UNC), and the 
Michigan State University (MSU). 
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of sight intersects three spiral arms (Sagittarius-Carina, Scutum- 
Crux, and Norma) . In addition to the tens of optical H ii re- 
gions identified by Rodgers et al. radio observations (e.g. 



ICaswell & Havnea 



19871) have revealed many other Hii regions 



heavily obscured by interstellar dust. The ionizing stellar clus- 
ters of som e of these objects have been studied in detail in the 
last d e cade jRoman-Lopes et alj|2003l: iRoman-Lopes & Abrahaii 



2004 iRusseil et alj |2005| ; Iskinner et al.l |2009|: Ichavam'a et al 



20101) . However, the ionizing stars of many of the Hii regions in 



the area have not been identified yet, and therefore their distances 
remain unknown. 

The 2:5 X 2:5 r adio Hn region GAL 331.31-00.34 
dCaswell & HavnesI Il987l) is one of the interesting and un- 
explored objects located in Norma. It is associated with the 
IRAS source 16085-5138 and with the 2:i x i:5 infrared 
cluster [DBS2003]157 centred at RA(J2000.0) = 16''12'"20% 
Dec.(J2000.0) = -51°46'14", which was identified by 



2 M. C. Pinheiro et al. 



iDutra et al] (|2003|) using data from the 2MAS S survey. Methanol 
(lEllingsen et alj Il996l : iKuchar & Cl^ Il99l IWalsh et alj [1993 : 
ICaswell et alJ bOOOh and hvdroxvl dCaswell et alj Il980h maser 
emiss ion indicates the existence of a massive star-forming re- 
gion ( IWalshetal.lll997l) . This is corrobo rated by the de t ection 
of CS(2-1) and SiO molecular lines ( [Bronfman et"aLl 1 19961: 
[Hariu et aL|[l998[ , respectively) which require H2 densities higher 
than 10" cm"^ 

GAL 33 1.3 1-00.34 is located at the east part of a large (1°) 
complex of Hn regions, which includes eight bright extended ra- 
dio sources ( [Amaral & Abrahanj [T991I) . A high-energy gamma- 
ray sourc e, HESS J1614-518, w as discovered at the edge of this 
complex ( [Aharonian et al.|[2006h . No counterpart has been found 
for this source among the most plausible classes of objects, like 
pulsar wind nebulae, supernova remnants. X-ray binaries, or ac- 
tive galactic nuclei. The superbubbles produced by the strong wind 
activity of OB associations could power high-energy gamma-ray 
sources jParizot et al.[ [20041) . Some other extended high-energy 
sour ces have also been fou nd associated with young stellar clus- 
ters | [ Aharonian et al. I [20071) . To figure out the energetics of these 
phenomena, it is essential to obtain the physical characteristics of 
the associated stellar population as well as an accurate estimate of 
its distance. 

In this work, we present JHK photometry of the stellar cluster 
[DBS2003] 157 and near-infrared (NIR) spectrophotometry of nine 
selected candidate stellar members. We aim to identify the ionizing 
sources of the associated Hn region GAL 331.31-00.34 and to es- 
timate its distance using the spectroscopic parallax method. In ad- 
dition to that, JHK photometry and mid-infrared (MIR) data from 
the GLIMPSE survey canied out with the Spitzer Space Telescope 
teeniamin et alj[2003h are used to identify young stellar objects in 
the area. 



2 OBSERVATIONS AND DATA REDUCTION 
2.1 Imaging photometry 

JHK photometric observations were performed at the Observatorio 
Pico dos Dias (OPD), Brazil, in April 2010. We used the NIR 
Camera CamlV attached to the 0.6-m Boiler & Chivens telescope 
to obtain frames with a field of view of 8' x 8' and spatial scale 
of 0'.'48pixer'. A Hawaii CCD detector of 1024 x 1024 pixels 
and a set of J, H, K^ filters were used. The K^, naiTow band fil- 
ter CI (FWHM ~ 0.023 /im), centred at 2.138^m, was chosen 
to avoid the contamination by B rackett-y nebular emission (see 
iRoman-Lopes & AbrahamI l2006a[) . To prevent saturation of the 
brightest stars in the J and H bands and to minimize the high 
thermal noise in the /fs-band, multiple short exposures were taken 
in each filter. The total integration times were 1260 s, 1575 s and 
5250 s in the H-, and Ci -bands, respectively. In order to subtract 
the background emission, images were obtained at five different po- 
sitions: the centre and its four adjacent positions, displaced between 
90" to 120". These large displacements were necessary to build a 
sky image uncontaminated by nebular emission. For each band, the 
final sky image represents the median value of each pixel taken 
over the five dithering positions. This technique allows the subtrac- 
tion of the background and at the same time preserves the nebular 
emission in the final images. The typical seeing ranged from 1" 1 
to r.'5. Dark exposures and dome flat-fields were also taken both at 
the beginning and at the end of each night. 

Standard procedures for NIR stellar CCD photometry in rel- 




Figure 1. CamlV //-band image of the cluster [DBS2003] 1 57. The ~ 2' x2' 
black square sun'ounds the cluster as indicated in DBS2()03. Blue marks 
denote all the sources found in the field by the daofind IRAF routine with 
a detection threshold of Act above the local background. The stars selected 
for the follow-up spectroscopy are numbered in red. 
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Figure 2. Top: Propagated photometric errors versus the /, //, and mag- 
nitudes. Bottom: 2MASS versus CamlV magnitudes. 



atively crowded fields were performed with the IRAlQ software. 
All frames were individually dark-subtracted and flattened, and 
the background contribution subtracted. Eventually, they were 
aligned and trimmed, to finally generate J-, H- and K-band im- 
ages by median-combining all frames of each filter. By choos- 
ing the median-filter, we diminished the influence of bad pixels 
and cosmic rays on the final image. This procedure resulted in 
i, H, and K^ "final" frames centred at RA(J2000) = 16''12'"24% 
Dec.(J2000) = -5r47'27", covering an area in the sky of ~5'x4'. 
Figure [T] shows one of these frames (//-band). The daofind IRAF 
routine was used in the determination of the physical coordinates 
of the stars in each frame, adopting an object detection threshold 
(in counts) of 4cr above the local background. 

The instrumental magnitudes were computed using the 
daophot IRAF routines. Since some parts of the cluster are too 
crowded for aperture photometry, the point spread functions (PSFs) 
were obtained using the p.s/ routine. Simultaneous PSF-fits for all 

' IRAF is distributed by the National Optical Astronomy Observatory, 
which is operated by the Association of Universities for Research in As- 
tronomy, Inc.. under contract to the National Science Foundation. 



the stars listed by daofind were performed with the allstar routine. 
A fitting radius of about 1 x {FWHM) and a PSF radius of ~ 3.5 
times this value was assumed for all stars. 

The photometric calibration was performed using data from 
the 2MASS survey. For this task, seventy-three no-blended stars 
in both 2MASS and CamlV frames were chosen outside the ~ 
2' X 2' rectangular region (see Fig.[D centred at the coordinates of 
[DBS2003] 157. Using the least squares method, we adjusted lin- 
ear transformation equations for converting the instrumental mag- 
nitudes to the 2MASS system. The standard deviations of these 
transformations, defined as 

'^' = ^J] ['"2MASS -Po-Pi m[iL.iv]' . (1) 

were 0.034, 0.026, and 0.036 for the J, H, and K, (CI) fil- 
ters, respectively. The final photometric errors were obtained as 
the quadratic sum of the transformation errors and of the indi- 
vidual error in the instrumental magnitudes given by the allstar 
task. No significant (J - H) and {H - K) colour-dependence be- 
tween 2MASS and CamlV photom etric systems has been found 
( iRoman-Lopes & Abraham|l2006lj) . Figure |2] shows the linear re- 
lation observed between the 2MASS and CamlV magnitudes and 
the dependence of the errors on the photometric magnitude for 
all stars detected. By analysing the magnitude histograms, log(A') 
vs. m, we estimate the completeness limits of the photometry as 
J = 15.5 mag, H = 15.2mag, and A^s = 13.5 mag, defined as 
the point where the log(A') vs. m plot deviates from a straight line. 
These values are similar to the 2MASS limits on the J and H bands 
(15.8 and 15.1 mag, respectively), but the A", completeness limit of 
the CamlV data is 0.8 mag lower. 

Table [ATI presents the results of the photometry. Columns 1 
lists our identification of the star; columns 2 and 3 provide the equa- 
torial coordinates obtained with the ccmap IRAF routine applied to 
pixel-coordinates listed by the daofind; columns 4, 5 and 6 list the 
H, and A^s photometric magnitudes and errors. 

At a spatial resolution twice as better as the 2MASS survey, 
we to identified many sources that were previously unresolved in 
the 2MASS data, which in turn improves significantly the photo- 
metric accuracy, especially for stars situated in the crowded areas. 

2.2 Spectroscopy 

NIR spectroscopic data of nine stars in the GAL 33 1.31- 
00.34 region were acquired with the Ohio State Infrared Im- 
ager/Spectrometer (OSIRISB coupled to the 4. 1-m telescope of the 
Southern Observatory for Astrophysical Research (SOAR), located 
at Cerro Pachon, Chilean Andes. Spectra in H-, and A'j-band 
were taken with OSIRIS in the low-resolution (R ^ 1200) multi- 
order cross-dispersed (XD) mode. In this mode, the instrument 
operates with a f/2.8 camera and a short (27"x 1") slit, covering 
the three bands simultaneously, in adjacent orders. The raw spectra 
were acquired using the standard AB nodding technique. Multiple 
short exposures were taken at each nod position, totalling 8-to-20 
individual frames, giving a "final" exposure time of 8-to-40 min- 
utes. Table [T] presents a journal of the spectroscopic observations. 
It lists our identification of the star, as defined in Table lAll the iden- 
tification in the 2MASS catalogue, the number and the duration of 
the individual exposures, and the observation dates. All the target 

^ OSIRIS is a collaborative project between the Ohio State University and 
Cerro Tololo Inter- American Observatory (CTIO). 
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Table 1. Journal of spectroscopic observations 





ID 


Indiv. exp. 


N° of 


Date 


This work 


2MASS 


time (s) 


exp. 




IRS 339 


J16122002-5146262 


120 


16 


2008/07/10 


IRS 298 


J16122053-5146460 


75 


12 


2011/05/08 


IRS 287 


J16122911-5146503 


75 


12 


2011/05/08 


IRS 497 


J16122925-5 147004 


120 


12 


2011/06/02 


IRS 355 


J16123324-5146173 


100 


08 


2011/06/02 


IRS 176 


J16122445-5 147492 


120 


12 


2011/06/03 


IRS 444 


J16122071-5147075 


90 


08 


2011/06/03 


IRS 289 


J16121561-5146504 


75 


20 


2011/06/12 


IRS 452 


J16122326-5146188 


120 


20 


2011/06/12 



stars were OB candidates selected from our photometry, with one 
exception: IRS 339, observed in 2008 and selected from an initial 
analysis based on the 2MASS data. 

A- and G-type spectroscopic standard stars were observed 
immediately before and after the "science" targets at similar air 
masses in order to remove telluric atmospheric absorption elfects 
from the "science" spectra. These raw spectra were reduced us- 
ing the CIRRED package and usual IRAF tasks. Two-dimensional 
frames were sky- subtracted for each pair of images taken at the 
two nod positions A and B, followed by dividing of the resultant 
image by a master flat. Thereafter, wavelength calibration was ap- 
plied using sky lines; the typical error (1 cr) for this calibration is 
estimate as ~ 12 A. The multiple exposures were combined and one- 
dimensional spectra were generated. 

Telluric atmospheric correction using the spectroscopic stan- 
dard stars completed the reduction process. In this last step, we 
divided each "science" spectrum by the spectrum of the AO V spec- 
troscopic standard star free of photospheric features, carefully re- 
moved by interpolating across their wings using the continuum 
points on both sides of the line. In the case of the //-band, the sub- 
traction of the hydrogen absorption lines could not be successfully 
done by this method because of the small separation between the 
multiple lines of the Brackett series and some strong telluric fea- 
tures. Therefore we proceeded to remove the hydroge n lines from 
the AO V //-band spectrum using the technique used bv lBlum et al.l 
( Il997h . Basically, one obtains a spectrum composed of only the 
profiles of hydrogen lines dividing the spectrum of a A-type stan- 
dard star by a spectrum of a G-type star, whose H i intrinsic lines 
(fairly weaker than in AO V line spectra) and other features of G- 
type spectra were previo usly removed by hand, using as template 
the NOAO solar atlas of iLivingston & Wallac3 ( Il99ll) . After that, 
these profiles are used to correct the AOV spectrum, generating 
an //-band spectrum free of the hydrogen Brackett lines. Finally, 
telluric bands were removed using the IRAF task telluric. This al- 
gorithm interactively minimizes the RMS in specific regions of the 
sample by adjusting the wavelength shifts and intensity scales be- 
tween the standard and "science" spectrum before performing the 
division. The wavelenght shifting corrects possible errors in the dis- 
persion zero-points, whereas the intensity scaling equalizes airmass 
differences and variations in the abundance of the telluric species. 
Typical values of the shifts varied around ~2 A, while the scaling 
factors were smaller than 10%. 
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Figure 3. Spectra of the hot stars identified in the field of GAL 331.31-00.34 (in black) and the comparison spectra for a 06 V star (HD 101190) and a 09 V 
star (HD 93028) (in red). 



by comparison with the 0.8-5 \im spectral library of cool stars by 
lRavneretal1 l l2009l) and are displayed in Fig. lAl [ (Appendix). They 
are probably field stars and will not be included in our discussion 
on the region. The other three spectra correspond to early-type stars 
(Figure [3}. 



The detection of Hen lines at 1.692^m and 2.188^m in the 
spectra of IRS 298 and IRS 339 restricts the spectral types of both 
of these stars to 09 or earlier. On the other hand, the relative 
strength of these lines compared with the He i features at 1 .279 yum, 
1.700 /jm and 2. 1 12 /jm indicates that IRS 339 is a late-0, whereas 
IRS 298 is an early-O star. The detection of the N iii line at 2. 1 15 yum 
in emission, together with the broad profile of the H i lines, indi- 
cates that these two objects are dwarf stars, even though we do not 
discard that IRS 339 may be a supergiant. The distinction between 
the various luminosity classes can be doubtful in the NIR because 
of the paucity of lines suitable for this purpose and the large uncer- 
tainties of their equivalent widths jHanson et alJI 199^120051) . 

Figure|4]presents a comparison between the equivalent widths 
(W)) of four selected He lines of dwarf stars, ranging from 06 to 
09.5 spectral types, along with the values measured for IRS 339 
and IRS 298. In this comparison, we included only lines with equiv- 
alent widths reasonably well determined and discarded Hi lines, 
whose profiles are highly dependent on the reduction process (see 
Sec. 12. 2t . Based on this scheme, the arguments presented above, 
and a visual comparison with other dwarf 0-star spectra obtained 
with the same instrumental configu ratio n and the s pectra of the at- 
lases bv lHanson et"!!] d 19961 . [200a) and lBik et al] i2005l) , we clas- 
sified IRS 298 as 06 V and IRS 339 as 09 V. 

The third stellar spectrum displayed in Fig. |3]is too noisy to 
be accurately classified. However, the strength of its He i lines and 
the tentative detection of Si, C, N, O suggest that IRS 335 is an 
09-to-B 1 type star. 
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Figure 4. Equivalent width W,i of He-lines versus the spectral type of dwarf 
stars (filled symbols) and a comparison with the values of W,i measured for 
IRS 298 and IRS 339 (open symbols, plotted at the spectral types 06 and 
09, respectively). 



3 RESULTS 

3.1 Spectral classiflcation 

He: optical a bsorption lines are well-known spectral features of 
OB stars (cf . : Walbom & Fitzpatrick 1990). If Hen absorption lines 
are also present, there is a substantial indicative of a n O-type star. 
Since such features are also found in NIR spectra jHanson et al.l 
Il996l) . the presence of these lines must still be the primary criterion 
to be observed when classifying hot star spectra in this band. Hi 
Brackett series could also be used. However, it would require spe- 
cial attent ion since these lines are highly contaminated by nebular 
emission ( iBik et al.ll2005l) . 

The spectra of six of the stars listed in Table [T] were found 
to correspond to late-type objects. They were classified visually 
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3.2 Spectroscopic parallax distance 

The distances of the classified stars IRS 298 and IRS 339 have been 
calculated from their apparent J and H magnitudes listed in Ta- 
ble |AT] using the method of spectroscopic parallax. The IR extinc- 
tion curve by[stead & Hoare (2009), which results in Ej-h/Eh-k = 
1.91, was adopted for the reddening corrections. As intrinsic stellar 
parameters, we used the absolute magnitude s for the various spec- 
tral types compilat ed by 



Vacca et al 



( 1 19961) and the non-reddened 



colour indices by iKoomnee i ( Il983l) . both transformed into the 
2MASS photometric system using the relations given bv lCarpenteJ 
( l200lh . According to this method, we obtained a heliocentric dis- 
tance ofd = 3.24 ± 0.56 kpc (Aj = 3.03, A v = 10.44) for IRS 339 
and d = 3.35 ± 0.61 kpc for IRS 298 {Aj = 3.51, Ay = 12.10). 
Based on the average distance of these two stars, we estimate 
the heliocentric distance for the Hii region GAL 331.31-00.34 of 
d = 3.29 ± 0.58 kpc and a visua l extinction of (Ay) ~ 1 1 mag. 

ICaswell & Havnej ( Il987l) obtained the radial velocity of 
Vlsr = -64kms"' for GAL 33 1.3 1-00.3 4, based on radio rec om- 
bination lines. According to the model o f iBrand & BUtl ( 1 19931) for 
the rotation curve of the Galaxy, this figure corresponds to the near 
and far kinematic distances of 4.2 kpc and 10.7 kpc, respectively. 
Thus this Hn region is definitely at the near kinematic distance. 



3.3 Inventory of the ionizing stars and the Lyman continuum 
luminosity 

In order to search other potential ionizing stars besides those 
confirmed spectroscopically, we reanalysed the NIR colour-colour 
(CC) and colour-magnitude (CM) diagrams (Figs.|5), assuming the 
distance of 3.29 kpc o btained in section 13.21 The intrinsic colours 
assumed are given by iKoomnee lill983h. while the absolute mag- 
nitudes in 7-band are from Vacca et al. ] d 19961) for 0-stars and 
from lWegnerl ( l2007h for the other spectral types, all of them trans- 
for med into th e 2MA SS photometric system by the relations given 
bv ICarpenteil (l200ll) . The classical T-Tauri star locus shown in 
the CC diagram is taken from iMever et"ar 1 ( Il997h . while the lo- 
cus of the Herbig Ae/Be c andidates is indicated in accordance 
with iLada & Adarni l ll992h . According to the CC diagram, we 
suspect that IRS 1 14 and IRS 446 are stars with spectral types be- 
tween BO V and B 1 V. However, they are heavily obscured by dust 
(Ay > 18 mag) and, in this region of the diagram, the classification 
is highly dependent on the slope adopted for the r eddening vector. 
iNishivama etal] l l2006l) and lStead & Hoard ( |2009|) showed that this 
slope depends on the direction of observation and, for large redden- 
ing, the scattering around the reddening vector might be very high. 
In these cases, a purely-photometric classification may not be very 
accurate. Thus, it could not be discarded that IRS 1 14 and IRS 446 
actually are pre-main sequence stars. 

An inventory of our search for hot stars in the field is pre- 
sented in Table 2, which lists the five stars identified as potential 
ionizing sources of GAL 33 1.3 1-00. 34 along with their photomet- 
ric and spectroscopic spectral classifications. 

To test the completeness of the set of ionizing stars found as- 
sociated with this Hii region, we compare the Lyman continuum 
photon luminosity A^Ly (photon s"') emitted by these stars with that 
needed to sup ply the radio contin uum flux. Adopting the Lyman 
luminosity bv lHanson et alj jl997h . these five stars would amount 
to a total of A'Ly = lO"*"" photons s"'. In fact, this total is strongly 
dominated by the contribution of the two 0-type stars. On the other 
hand, a Lyman continuum luminosity of A^Ly = lO'*^ ''' photons s"' 
is inferred from the radio continuum flux density Sy at 22 GHz 
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potential ionizing stai's are also indicated. 



Table 2. Ionizing stars identified in GAL331. 31-00.34 





Star ID 


Classifications 


This work 


2MASS 


Photometrical 


Spectral 


IRS 298 


J16122053-5146460 


mid-O 


06 V 


IRS 339 


J16122002-5146262 


late-O 


09V 


IRS 355 


J16123324-5146173 


BO: 


BO: 


IRS 1 14 


J16123539-5148304 


Bl: 




IRS 446 


J16120863-5 146482 


BO: 
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Table 3. Young Stellar Objects candidates from the NIR photometry 
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Figure 6. MIR colour-colour diagram of star s in a region with 2^5 of radius 
and centre matching with //-band frame presented in Fig. [T] YSOs with 
MIR colour excesses he rightward of the dotted line. Class III YSOs and 
back/foreground objects must fall in the clump around the (0, 0) point. 



measured bylAmaral & Abrahaml(ll99lh . using the expression: (see 
lRubiij|l968[) 



5- 59x10^^ / V \°-' 0.45, .2 

e+/(H+ -)-He+)) \5GHz/ 



(2) 



1 -h,/;<He+/(H+ -)-He+)> 

where we assume the abundance ratio of He^/H+ = 0.5, a typical 
electron temperature = 8000 K, and / = 0.65 for the fraction of 
He-recombination photons energetic en ough to ionize the H . Using 
the flux densities at 5 GHz m easured bv lCaswell & HavnesI l ll987l) 
and iKuchar & Cla3 i ll 9971) almost identical results are found. 
Thus, we conclude that the O-type stars IRS 298 and IRS 339 can 
cope with the ionization of this nebula. The three stars photomet- 
rically classified as BO-to-B 1 V type and other cooler stars would 
not contribute significantly to the ionization of the nebula. 

3.4 YSO candidates 

The intrinsic IR excess found in some stellar objects is an evidence 
of accreting disks or iirfalling envelopes. The identification of a 
group of objects with these properties indicates the existence of in- 
tense and recent star formation, which makes these regions interest- 
ing places to search for pre-main-sequence stellar objects (PMS), 
for example. The analysis of NIR and MIR colour-colour diagrams 
allows us to identify PMS stars and infer their natu re based on 
their p hotometric properties. For example, according to lAllen et al.l 
( l2004h , objects with accreting disks, named Class II YSO, occupy 
the locus of reddened classical T-Tauri in NIR colour-colour dia- 
grams, whereas objects whose emission is dominated by infalling 
envelopes, named Class I, have more reddened colours. YSOs with 
weak or no IR-excess are named Class III, and are not discriminated 
in CC diagrams, since they are found along with ZAMS stars. 

The analysis of the NIR colour-colour diagram presented in 
Fig. |5] indicates a total of 36 YSO candidates, which are listed in 
Table [3] together with their tentative classification between T-Tauri 
(TTau) or Herbig Ae/Be (HAeBe) stars. In this list, IRS 446 and 
IRS 114 could also be main sequence stars very obscured by dust. 
IRS 150 is the most reddened object detected in this work, with 
Ay > 20 mag. Considering the error bars of the photometry and 
extinction laws, IRS 150 has colours of a very reddened T-Tauri 
star. Redwards in the CC NIR diagram, IRS 478 presents typical 





Identifier 






L.iass. 


This work 


2MASS 


GLIMPSE 


IRS 42 


J16123385-5149269 


G331 


3195-00.3979 


HAeBe 


IRS 94 


J16 122500-5 148459 






HAeBe 


IRS 114 


J16123539-5148304 


G331 


3330-00.3892 


TTau+ 


IRS 133 


J16121385-5148204 


G331 


2944-00.3491 


HAeBe 


IRS 138 


J16 122037-5 148 157 


G331 


3076-00.3597 


HAeBe 


IRS 150 


J16 123309-5 148076 


G331 


3331-00.3806 


TTau 


IRS 172 


J16120964-5147550 






TTau 


IRS 178 


J 16123437-5 147487 


G331 


3391-00.3790 


HAeBe 


IRS 184 


J16121891-5147434 


G331 


3110-00.3506 


HAeBe 


IRS 197 


J16121395-5147382 


G331 


3025-00.3409 


HAeBe 


IRS 224 


J16121561-5147255 






HAeBe 


IRS 246 


J16121386-5147139 


G331 


3071-00.3357 


HAeBe 


IRS 256 


J16 122966-5 147088 






HAeBe 


IRS 283 


J16123053-5146525 






HAeBe 


IRS 330 


J 16122692-5 1463 14 






HAeBe 


IRS 337 








HAeBe 


IRS 344 


J16122121-5146255 






HAeBe 


IRS 362 


J16122064-5 146136 






HAeBe 


IRS 376 


J16121207-5 146026 






HAeBe 


IRS 383 


J16121868-5145597 






HAeBe 


IRS 390 


J16123848-5145538 


G331 


3687-00.3630 


TTau 


IRS 402 


J16 123964-5 145474 






HAeBe 


IRS 445 


J 16122483-5 146507 






HAeBe 


IRS 446 


J16120863-5 146482 


G331 


3022-00.3213 


TTaut 


IRS 454 


J16122704-5146168 


G331 


3427-00.3474 


TTau 


IRS 478 


J16 122996-5 148 154 


G331 


3258-00.3766 


HAeBe 


IRS 481 


J16122752-5148042 


G331 


3232-00.3700 


TTau 


IRS 483 


J16121153-5147565 


G331 


2947-00.3402 


HAeBe 


IRS 484 


J16122607-5 147546 


G331 


3223-00.3655 


HAeBe 


IRS 487 


J1612261 1-5147500 


G331 


3232-00.3646 


TTau 


IRS 488 


J16120772-5147483 


G331 


2891-00.3319 


HAeBe 


IRS 5 10 


J161221 12-5146395 


G331 


3274-00.3416 


TTau 


1RS511 


J16120732-5146382 






HAeBe 


IRS 524 


J16123883-5148207 


G331 


3413-00.3932 


HAeBe 


IRS 531 


J16120752-5146336 


G331 


3028-00.3164 


HAeBe 


IRS 533 








HAeBe 



fSee Section[33] 

Table 4. Identifiers of the YSOs candidates from the MIR GLIMPSE pho- 
tometry 



This work 



Surveys 



MIR 



NIR 



GLIMPSE 



2MASS 



CEl 




G33 1.2999-00.3784 


J16122314-5149240 


CE2 




G33 1.3093-00.3762 


J16122521-5148551 


CE3 




G331.3205-00.3817 




CE4 


IRS 552 


G331.3131-00.3653 


116122341-5148167 


CE5 


IRS 168 


G33 1.3298-00.3735 


116123029-5147575 


CE6 


IRS 176 


G331.3203-00.3615 


116122445-5147492 


CE7 




G331.3341-00.3619 


J16122844-5147163 


CE8 


IRS 645 


G331.3541-00.3664 


116123527-5146387 


CE9 


IRS 355 


G33 1.3543-00.3585 


116123324-5146173 


CEIO 




G331.3416-00.3464 


J16122648-5146177 


CEll 




G33 1.3497-00.3497 
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Figure?. Spectral energy distributions (between 1.25 and 10//m) of objects 
with MIR colour excess. Photometric data extracted from the 2MASS and 
GLIMPSE surveys. 
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Figure 8. Spectral energy distiibution (between 1.25 and 100 /jm) of the 
source CE 10 based on data extracted from the 2MASS, GLIMPSE, and 
IRAS surveys. The dotted line (with spectral index of a » 3.6) represents 
the fitting of the SED between 2 and 25 )im. 



1.0. However, from the MIR photometry alone, we are not able to 
decide whether or not these objects are low-mass YSOs. 

Figure [7] presents the spectral energy distributions (SEDs) of 
the MIR sources classified as YSOs, obtained by combining their 
NIR and MIR magnitudes. Sources CE 1 and CE 10 present a high 
spectral index a = dlog(AFA)/dlogA near A ~ 2yum, wh ich is 
chara cteristic of protostars in an early stage of evolution jLadal 
Il987h . 



3.5 The nature of IRAS 16085-5138 

IRAS 16085-5138 is the only IR AS source in the field Acc ord- 
ing to the classification scheme of IWood & Churchwelll ( Il989l) . its 
IRAS colours are compatible with an ultra-compact Hit region. 
Methan ol and hydroxyl ma s ers ha ve been detected at its vicin- 
ity (e.g. ICaswell etaLlll980l . l2000h . We identify this source with 
CE 10, which presents the highest spectral index among all the CE- 
objects identified in this work, and is the only one detected inside 
the IRAS position error ellipse. Using photometric data from the 
2MASS, GLIMPSE, and IRAS surveys, w e obta i ned th e 2-100 yum 
SED of this source (Fig.[8]l. According to lLadal l ll987h . a spectral 
energy distribution rising longward of A = 2 yum, such as that shown 
by this object, is characteristic of Class I objects. We evaluated the 
spectral index between 2 and 25 (the dotted line in Fig. (8)1, and 
found a » 3.6, typical of YSO immersed in a dense protostellar 
envelope, which is also in accordance with the classification drawn 
from the MIR colour-colour diagram. 

Since most of the energy of this object is emitted in the in- 
frared, the bolometric luminosity of the embedded protostar can 
be approximated by integrating the SED between 2 and lOOyum. 
Adopting the distance of 3.29 kpc and Ay = 10 mag, we obtain a 
lumi nosity of L = 7.7 x IQ-^Lq. Since pre-main-sequence models 
(e.g. llbenlll965l) indicate that early-type stars evolve at constant lu- 
minosity from the Hayashi track up to the main sequence, this result 
implies that the lower limit for this stellar mass is M > 10 Mp. This 
mass value con'esponds to a BO-B 1 ZAMS star (cf. iHanson et al.l 
Il997h . 



colours of a very embedded YSO in an early stage of formation. 
A third group of stars appears below the T-Tauri locus, showing 
(H - Ks) > 0.4 mag. In this region, Herbig Ae/Be stars, which are 
PMS star s more massive than T -Tauri, can be confused with Class 
I sources ( lLada& Adamslll992h . galaxies or other background ob- 
jects, making a classification based only on photometric data unre- 
liable. 

Figure |6] presents a MIR colour-colour diagram made with 
data from the GLIMPSE survey. This diagram is used to identify 
and classify objects with MIR colour excess (hereafter denoted as 
CE). Within a region of 2^5 radius centred at the NIR //-frame 
(see Fig.[T), a total of 58 sources were simultaneously detected in 
the four GLIMPSE/IRAC-bands. We identify 1 1 MIR sources with 
colour excesses, which are listed in Table |4] Some of them are so 
obscured by dust that they have not been detected in the NIR bands. 

Low-mass objects with M IR colour excess can be classified 
according to scheme by I Allen e t al. (2004). Class III objects, along 
with back/foreground objects, have the MIR colours [3.6] - [4.5] ~ 
and [5.8] - [8.0] ~ 0, Class II have colours within < [3.6] - 
[4.5] < 0.8 and 0.4 < [5.8] - [8.0] < 1.1, whereas Class I objects 
are expected to present colours [3.6]-[4.5] > 0.7 and [5.8]-[8.0] > 



3.6 MIR shell structure 

Figure |9] presents a RGB combination of 8.0;um (R), 5.8yum (G), 
and 3.6yum (B) GLIMPSE images that show a broken shell-like 
structure with a diameter of » 4' associated with GAL331.31- 
00.34. At the distance of 3.29 kpc, this shell would have a diameter 
of a; 3.8 parsec. This shell, identified as [CPA200 6] S62 in the SIM- 
BAD d atabase, is named S 62 in the catalogue by I Church well et akl 
l|2006f). which contains over 300 MIR arcs. The authors interpreted 
the arcs as the projection on the sky of three-dimensional bubbles, 
and found that only one quarter of them would be associated with 
H n regions. In these cases, the bubbles would be expanding driven 
by the stellar wind and the radiation pressure from young massive 
stars. These stars would also power the bubbles, which are traced 
by the polycyclic aromatic hydrocarbon (PAH) emission produced 
in the photodissocia tion regions located at t he edge of the Hii re- 
gions. According to iDeharveng et al.l ( l2010h . most of the bubbles 
are in fact associated with H ii regions, and more than a quarter of 
them could be triggering the formatio n of a new generation o f stars, 
as seems to be the case of RCW 120 ( IDeharveng etai]|2009l) . 

In the case of GAL 331.31-00.34, the IRAS source 16085- 



iKuchar & Cl^ 1 19971; IWalsh et alj Il998l ; ICaswell et all 



dEllinesen et al. 


1996 


Caswell et al.l 


200C 
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Figure 9. RGB combination of 8.0^(m (R), 5.8^im (G), and 3.6^im (B) 
GLIMPSE images. The red circles indicate tlie objects showing MIR colour 
excess hsted in Table|4] In these three MIR bands, the peak emission coin- 
cides with the source CE 10, which we identified as the GLIMPSE counter- 
part of IRAS 16085-5138. 

119801) . typical tracers of ongoing massive star formation, are lo- 
cated at the bright northern border of the shell. At these posi- 
tions, CS(2-1) and SiO molecular emissio n, which require densi- 
ties higher than 10^_cm"^, were detected bv lBronfman et aL l( ll996l) 
and lHariu et al.l ( Il998h , respectively, who carried out radio surveys 
towards IRAS sources with colours of UC-HII regions. This possi- 
bly indicates that the expansion of the northern part of the bubble 
has been hampered by a dense molecular environment. The same 
phenomenon may also occur in other parts of the bubble, but addi- 
tional observations of CS and NH3, which are good tracers of dense 
molecular regions, are still lacking. 

We also find that the T Tauri and the Herbig Ae/Be candidates 
found from the NIR photometry are found preferentially towards 
the bubble, while the YSO candidates found from the MIR data are 
located in a limited north-south band crossing the openings of the 
broken shell. 




le'ia^ae' 30' 24' is' 12" 

a (J2000) 



Figure 10. Contour levels of the IRAC 8.0/jm image superimposed on H- 
band image. The yellow triangles show the positions of the methanol and 
hydroxyl maser sources. The cyan and orange triangles locate the low-mass 
and high-mass YSO candidates (possible reddened T-Tauri and the Herbig 
Ae/Be stars, respectively) classified according to the NIR photometry and 
listed in Table[3] The green circles indicate the objects showing MIR colour 
excess listed in Table |4] The blue dots point to the hot stars spectroscopi- 
cally identified. 

and 27 Herbig Ae/Be candidates. From the GLIMPSE data, 1 1 ob- 
jects with MIR colour excesses have been found. 

iii. The MIR counterpart of the IRAS source 16085-5138 was 
identified and its spectral energy distribution between 2 and 100 /jm 
was analysed. We concluded that IRAS 16085-5138 is a massive 
YSO in an early stage of formation, with luminosity L ^ 7.7 x 

10' Lq and 

mass M > W Mq. 

High spatial resolution observations in radio frequency of 
the ionized and molecular gas are necessary for deeper studies 
of the dynamics and evolution of star formation in the field of 
GAL331. 31-00.34. Likewise, the possible connection between the 
diffuse high energy gamma-ray emission HESS 11614-518 and 
powerful winds and radiation fields of OB stars in the region should 
be subject of future hydrodynamic studies. 



4 SUMMARY 

We carried out JHK photometry in the direction of the Galactic 
Hii region GAL 331.31-00.34. We classified spectroscopically the 
main ionizing sources and estimated their distance. As a result, we 
identified other potential ionizing stars and performed their respec- 
tive photometric classifications. Data from the 2MASS, GLIMPSE, 
and IRAS surveys were also explored. Our main findings are: 

i. We identified, classified spectroscopically, and estimated the 
spectroscopic parallax distances of two O-type stars associated 
with this Hi: region: IRS 298 (06V, d = 3.35 ± 0.61 kpc) and 
IRS 339 (09 V, d = 3.24±0.56 kpc). Adopting the average distance 
of 3.29 ± 0.58 kpc and comparing the Lyman continuum luminosity 
of these stars with that required to ionize the nebula, obtained from 
radio continuum observations, we concluded that these two stars 
are the ionizing sources of GAL 331.31-00.34. 

ii. By analysing the NIR colour-colour diagram, 36 pre-main se- 
quence (PMS) objects could be identified and classified: 9 T-Tauri 
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Table Al. JHK photometry 



IRS g(J2QQQ) a'(J2QQQ) J H 



3 


16 


12:33.028 


-51:49:23.11 


16.166+0.070 


13.485±0.018 


12.245±0.065 


41 


16 


12:26.165 


-51:49:27.27 


13.479±0.028 


11.951±0.016 


11.199±0.037 


42 


16 


12:33.848 


-51:49:26.83 


13.792±0.031 


13.405±0.014 


13.045±0.1 12 


79 


16 


12:14.127 


-51:49:03.79 


1 1.658+0.029 


9.809 ±0.015 


9.021 ±0.015 


87 


16 


12:37.335 


-51:48:57.01 


13.042+0.025 


1 1.841±0.012 


1 1 .447 ±0.039 


94 


16 


12:24.999 


-51:48:45.83 


14.955±0.031 


14.212±0.014 


13.090±0.231 


98 


16 


12:39.840 


-51:48:41.18 


9.475 ±0.032 


8.185 ±0.013 


7.682 ±0.014 


102 


16 


12:23.151 


-51:48:40.40 


9.738 ±0.028 


9.044 ±0.008 


8.824 ±0.015 


108 


16 


12:34.752 


-51:48:33.93 


15.580+0.039 


13.894±0.017 


12.969±0.085 


113 


16 


12:31.255 


-51:48:30.52 


12.930±0.028 


11.270±0.009 


10.501 ±0.0 18 


114 


16 


12:35.383 


-51:48:30.43 


15.455+0.038 


13.179±0.018 


11.801±0.045 


117 


16 


12:28.521 


-51:48:29.84 


12.435±0.029 


11.458±0.014 


11.163±0.022 


133 


16 


12:13.855 


-51:48:20.25 


14.057±0.028 


13.593±0.022 


13.071±0.112 


135 


16 


12:38.748 


-51:48:16.20 


10.923±0.029 


10.217±0.013 


9.954 ±0.014 


138 


16 


12:20.382 


-51:48:15.79 


14.586+0.025 


14.140±0.024 


13. 745 ±0.195 


150 


16 


12:33.097 


-51:48:07.57 


15.624+0.053 


12.155±0.015 


9.842 ±0.013 


164 


16 


12:13.044 


-51:47:59.20 


13.887±0.029 


12.279±0.016 


11.543±0.028 


168 


16 


12:30.291 


-51:47:57.57 


1 1.302+0.026 


10.720±0.015 


10.448 ±0.0 13 


169 


16 


12:31.902 


-51:47:57.47 


10.657+0.028 


9.965 ±0.013 


9.775 ±0.016 


172 


16 


12:09.542 


-51:47:54.88 


15.896+0.050 


14.828±0.042 


13.973±0.360 


176 


16 


12:24.451 


-51:47:49.26 


14.321±0.026 


12.807±0.015 


12.347 ±0.081 


178 


16 


12:34.365 


-51:47:48.51 


13.576+0.029 


13.295±0.010 


12.752±0.131 


184 


16 


12:18.908 


-51:47:43.46 


14.291 ±0.026 


13.787±0.014 


13.393±0.197 


195 


16 


12:40.537 


-51:47:37.13 


1 1.914±0.032 


10.653±0.012 


10.148 ±0.015 


197 


16 


12:13.960 


-51:47:38.10 


13.89O±0.O3O 


13.448±0.024 


12.800±0.1 15 


205 


16 


12:36.513 


-51:47:32.56 


12.382±0.027 


1 1.490±0.016 


11.162±0.021 


217 


16 


12:19.360 


-51:47:28.72 


9.105 ±0.025 


8.155 ±0.011 


7.927 ±0.014 


224 


16 


12:15.507 


-51:47:25,57 


15.415±0.037 


14.917±0.029 


13.909±0.273 


230 


16 


12:29.728 


-51:47:19.59 


14.678±0.036 


11.831±0.016 


10.561 ±0.01 8 


234 


16 


12:12.487 


-51:47:18.36 


1 1.570±0.023 


10.803±0.006 


10.559±0.023 


241 


16 


12:36.560 


-51:47:15.40 


15.81 1±0.051 


13.884±0.022 


12.785±0.115 


246 


16 


12:13.859 


-51:47:14.03 


15.552±0.033 


14.650±0.016 


13. 543 ±0.275 


255 


16 


12:11.494 


-51:47:09.36 


15.456±0.080 


14.015±0.028 


13. 145 ±0.090 


256 


16 


12:29.559 


-51:47:08.77 


15.042±0.057 


14.215±0.043 


12.71 0±0.094 


258 


16 


12:37.134 


-51:47:07.52 


1 1 .005±0.029 


10.247±0.009 


9.977 ±0.019 


262 


16 


12:32.924 


-51:47:05,09 


12.240±0.026 


11.318±0.010 


10.978±0.021 


264 


16 


12:28.519 


-51:47:04.38 


11.391±0.026 


11.089±0.012 


1 1.095±0.025 


268 


16 


12:30.324 


-51:47:01.98 


11.528±0.024 


10.457±0.014 


10.035±0.015 


277 


16 


12:26.742 


-51:46:57.38 


1 1.536±0.026 


11.224±0.016 


11.200±0.031 


278 


16 


12:34.500 


-51:46:56.74 


10.497 ±0.027 


9.805 ±0.013 


9.514 ±0.017 


283 


16 


12:30.538 


-51:46:52.58 


15.152±0.085 


15.411±0.050 


14.345±0.31 1 


285 


16 


12:27.735 


-51:46:50.82 


15.592±0.056 


13.071±0.010 


12.045±0.043 


286 


16 


12:19.525 


-51:46:51,10 


15.052±0.025 


12.684±0.015 


11.585±0.041 


287 


16 


12:29.115 


-51:46:50.29 


1 1.203±0.025 


10.159±0.009 


9.736 ±0.012 


289 


16 


12:15.521 


-51:46:50.44 


14.330±0.021 


11.136±0.012 


9.643 ±0.013 


295 


16 


12:32.209 


-51:46:46.38 


15.284±0.037 


1 1.904±0.01 1 


10.302±0.013 


298 


16 


12:20.539 


-51:46:46.07 


1 1 .674±0.024 


10.288±0.O08 


9.593 ±0.018 


319 


16 


12:27.952 


-51:46:35.74 


10.950±0.024 


9.707 ±0.012 


9.216 ±0.013 


323 


16 


12:40.287 


-51:46:34.17 


14.192±0.026 


12.371±0.017 


11.540±0.034 


324 


16 


12:09.102 


-51:46:35.02 


1 1.319±0.024 


10.520±0.009 


10.205±0.013 


330 


16 


12:26.925 


-51:46:31.44 


14.554±0.025 


13.945±0.017 


13.499±0.192 


337 


16 


12:19.501 


-51:46:27.59 


14.956±0.078 


14.479±0.083 


13.584±0.215 


339 


16 


12:20.031 


-51:46:26.20 


1 1 .784±0.024 


10.596±0.007 


10.001±0.015 


344 


16 


12:21.212 


-51:46:25.36 


15.906±0.059 


15.476±0.055 


14.138±0.355 


355 


16 


12:33.242 


-51:46:17.41 


1 1.728±0.022 


10.691 ±0.015 


10.184±0.015 


362 


16 


12:20.553 


-51:46:13.75 


14.693±0.033 


14.312±0.023 


13.494±0.192 


363 


16 


12:32.580 


-51:46:13.00 


12.597 ±0.026 


11.052±0.014 


10.408 ±0.0 18 


369 


16 


12:36.212 


-51:46:08,58 


11.911±0.029 


10.73O±0.O13 


10.264±0.018 


374 


16 


12:37.823 


-51:46:03.37 


14.042±0.026 


12.130±0.012 


11.320±0.031 


375 


16 


12:35.031 


-51:46:02.14 


1 1.553±0.022 


10.814±0.014 


10.575±0.015 


375 


16 


12:12.069 


-51:46:02.71 


14.307±0.026 


14.016±0.018 


13.579±0.240 


383 


16 


12:18.702 


-51:45:59.54 


15.204±0.042 


14.690±0.032 


14.178±0.361 


390 


16 


12:38.478 


-51:45:53.96 


13.809±0.O27 


12.774±0.017 


12.078±0.047 


391 


16 


12:27.182 


-51:45:53.34 


13.145±0.023 


11.842±0.014 


11.218±0.027 


402 


16 


12:39.547 


-51:45:47.38 


15.756±0.057 


15.157±0.039 


14.322±0.353 


410 


16 


12:35.158 


-51:45:44.50 


15.733±0.046 


12.333±0.009 


10.799±0.022 


415 


16 


12:32.059 


-51:45:39,79 


1 1.884±0.025 


11.189±0.014 


10.925±0.020 


415 


16 


12:23.793 


-51:45:40.10 


1 1.733±0.023 


1 1.057±0.012 


10.795±0.023 


425 


16 


12:10.705 


-51:48:53.53 


13.672±0.026 


12.021 ±0.01 3 


11.345±0.031 






12'07 291 


-51'47'22 36 


12 806±0 019 


1 1 696+0 008 


1 1 330±O 017 


444 


16 


12:20.717 


-51:47:07.56 


11.615±0.035 


10.535±0.052 


10.144±0.050 


445 


16 


12:24.838 


-51:46:50.52 


14.336±0.032 


13.553±0.053 


12.820±0.143 


445 


16 


12:08.534 


-51:46:48.29 


15.899±0.085 


13.180±0.019 


11.565±0.033 


452 


16 


12:23.287 


-51:46:18.82 


14.254±0.045 


11.793±0.016 


10.557±0.017 


454 


16 


12:27.045 


-51:46:16.80 


14.315±0.033 


13.110±0.019 


12. 295 ±0.065 


478 


16 


12:29.948 


-51:48:15.55 


15.81O±0.10O 


13.951 ±0.075 


12.172±0.055 


481 


16 


12:27.537 


-51:48:04.28 


15.985±0.101 


14.350±0.055 


13.074±0.134 


483 


16 


12:11.535 


-51:47:56.37 


13.732±0.024 


13.250±0.014 


12.797±0.098 


484 


16 


12:26.065 


-51:47:54.50 


15.870±0.058 


15.093±0.097 


13.571 ±0.221 


487 


16 


12:26.122 


-51:47:50.07 


15.091 ±0.03 3 


13.783±0.022 


12.938±0.117 


488 


16 


12:07.740 


-51:47:48,32 


14.956±0.058 


14.385±0.053 


13.508±0.199 


493 


16 


12:28.865 


-51:47:21.05 


15.241 ±0.1 52 


14.295±0.154 


13.235±0.137 


494 


16 


12:12.964 


-51:47:13,43 


12.650±0.023 


11.589±0.011 


11.183±0.021 


497 


16 


12:29.250 


-51:47:00.40 


13.407±0.047 


11.917±0.043 


11. 351 ±0.052 


510 


16 


12:21.101 


-51:46:39.48 


15.272±0.150 


14.472±0.032 


13.269±0.110 


511 


16 


12:07.325 


-51:46:38.31 


15.433±0.053 


14.657±0.054 


13.508±0.143 


515 


16 


12:19.293 


-51:46:07.20 


12.907±0.028 


11.880±0.022 


11.509±0.035 


524 


16 


12:38.811 


-51:48:20.72 


13.954±0.139 


13.450±0.127 


12.917±0.109 


530 


16 


12:27.472 


-51:46:54.54 


15.507±0.153 


13.005±0.045 


12.004±0.061 


531 


16 


12:07.529 


-51:46:33.57 


14.01 1±0.018 


13.496±0.019 


13.037±0.135 


533 


16 


12:22.212 


-51:46:10.27 


15.790±0.058 


15.156±0.026 


14.192+0.242 



APPENDIX A: 



The young stellar cluster DBS2003 157 1 1 




\(//m) \{fj,m.) \{fj,m) 



Figure Al, 



. Cold stars in the field of GAL331. 31-00.34. 



